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We present results of a calculation of the electromagnetic pion form 
factor within the framework of QCD Sum Rules with nonlocal condensates, 
using a pcrturbative spectral density which includes 0(a s ) contributions. 



1. Introduction 

The pion — its distribution amplitude and form factor — stands tall as a 
role model for the modern description of hadrons in terms of quarks and 
gluons within QCD. At high momenta Q 2 , the pion form factor can be 
written as a convolution F 1T {Q 2 ) = ip° ut ® T(Q 2 ) ® ip™ on account of the 
factorization theorem, where the symbol (g) means an integration over the 
longitudinal momenta of the quark and antiquark in the pion state factor- 
ized at some scale [i 2 . All binding effects due to the nonperturbative color 
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dynamics are absorbed into the pion distribution amplitudes c^™ for the in- 
coming and </j° ut for the outgoing pion. However, at low and intermediate 
momenta Q 2 this factorization procedure becomes inapplicable because the 
long-range interactions cannot be separated into factorizing pion distribu- 
tion amplitudes ip-^, as above, so that the convolution approach is not very 
reliable. 

In a recent paper [I], we have proposed a different theoretical framework 
for the calculation of the pion form factor, which is based on QCD sum 
rules with nonlocal condensates (NLC) [21 El [4], employing an Axial- Axial- 
Vector (AAV) correlator. Suffice it to say here that our scheme differs 
from previous ones in several respects that will be considered in the next 
section. A prolegomenon of the results we derived (see Sec. [3]): the proposed 
method yields predictions for the spacelike pion form factor that compare 
well with the data in that momentum region which is currently accessible 
to measurements, covering also the range of momenta to be probed by 
the 12 GeV upgraded CEBAF accelerator at JLab. Our conclusions are 
summarized in Sec. HI 

2. Three-point QCD sum rule for the pion form factor 

Let us make here a second, more detailed pass, on some of the topics 
briefly mentioned in the Introduction, recalling the AAV correlator 



where q denotes the photon momentum (q 2 = — Q 2 ) and p 2 is the outgoing 
pion momentum. The quantities J M (x) = e u u(x)^^u(x) + edd{x)^/ ti d{x) and 
Jbu{x) = d(x)'j^ a u(x) , J^(x) = u{x)^ r ypd{x) are the electromagnetic 
current and the axial- vector currents, respectively, where e u = 2/3 and 
erf = —1/3 stand for the electric charges of the u and the d quarks. Referring 
for further details to [5l [6] for the case of local condensates and to [I] for 
nonlocal ones, we proceed by writing down the sum rule we will employ: 



//< 



d 4 xd 4 ye^ x -^{0\T J+Jy)J»(x)J 5a (0) |0> 



(1) 



d Sl ds 2 p 3 ( Sl ,s 2 ,Q 2 )e-^ +s ^ M2 







+<MQ 2 ,M 2 ) + cD (g - g) (Q 2 ,M 2 ). 



(2) 



Note that the quark-condensate contribution 



<J> (9 - 9) (Q 2 ,M 2 ) = $ 4Q (Q 2 ,M 2 ) + $ 2V (Q 2 ,M 2 ) + $ g - Ag (Q 2 ,M 2 ) (3) 
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contains the four-quark condensate (4Q), the bilocal vector-quark conden- 
sate (2V), and the antiquark-gluon-quark condensate (qAq), while the term 
*&g{Q 2 -,M 2 ) represents the gluon-condensate contribution to the sum rule. 
The crucial quantity in the sum rule is the three-point spectral density 



pf\si,s 2 ,Q 2 } 



Pl'{S\,S2,Q )+ ^ Ap y 3 J {si,S 2 ,Q ) 



(4) 



The leading-order spectral density has been calculated long ago [SJ [6], 

whereas the analogous next-to-leading order (NLO) version Apg (si, s 2 , Q 2 ) 
has been derived recently in [7\. The contribution from higher resonances 
is usually taken into account in the form 

Phr(si, s 2 ) = [1 - 0(s 1 < s o )0{s 2 < s )\ p 3 (si,s 2 , Q 2 ) (5) 

and contains the continuum threshold parameter sq. In the investigation [I], 
reported upon here, we use in the perturbative spectral density a version of 
the running coupling that avoids Landau singularities by construction (see 
for reviews in [HO [TO]). At the one- loop level, one has [8] 

asiQ2) = ? (mqvWd) " q 2 -^) (6) 

with b = 9 and A QC d = 300 MeV. 

The key elements of our analysis are these: (i) All interquark distances 
in the quark-gluon-antiquark condensate are nonlocal and the nonlocality 
is parameterized via the quark-virtuality parameter A 2 , [2] with the value 
0.4 GeV 2 . (ii) A modified Gaussian model for the nonlocal condensate is 
used, recently proposed in pj~j, and the prediction for is compared with 
the results derived in [12] , obtained by using the minimal Gaussian model [3l 
[T3l [T3] , and from other theoretical models |151 [TO] as well. The virtue of the 
modified NLC model is that it helps minimizing the transversality violation 
of the two-point correlator of vector currents, (iii) A spectral density is 
used that includes terms of 0(a s ), i.e., NLO perturbative contributions. 
Moreover, the coupling entering the spectral density is analytic, as stated 
above, so that the calculation of the pion form factor is not influenced by 
the Landau pole. It was shown in [TTJ, [12] that the Landau pole can obscure 
the predictions for the pion form factor even at momentum values much 
larger than Aq C d- 



3. Predictions 



Using the presented scheme, we obtain predictions for the pion form 
factor [1] , which are shown in Fig. [TJ Our results are presented in the form 
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Fig. 1. Scaled pion form factor Q 2 F n (Q 2 ) for the minimal (dashed blue lines) and the 
improved (solid blue lines) NLC model using A 2 = 0.4 GeV 2 in comparison with experi- 
mental data of the Cornell [18] (triangles) and the JLab Collaboration [19] (diamonds). 
The shaded band within dashed lines shows the minimal model, whereas the variance of 
the improved model is indicated by solid lines. The lattice result of [20] is also shown as 
a dark-grey strip. The two broken vertical lines mark the region, where the influence of 
the particular Gaussian NLC model used is not severe. 

of shaded bands, the aim being to include the inherent theoretical uncer- 
tainties of the method. The band contained within the solid lines gives the 
predictions we obtained with the improved Gaussian NLC model, whereas 
the corresponding findings from the minimal NLC model are shown in the 
band with the dashed boundaries. The central curves in each band (illus- 
trated by a thick solid and a thick dashed line) are interpolation formulas, 
whose explicit form can be found in pQ. For the understanding of these 
predictions it is instructive to remark that our method provides results that 
remain valid — though with reduced accuracy — even at higher values of the 
momentum up to Q 2 ~ 10 GeV 2 . 

Let us now discuss some technical details. The value of sq{Q 2 ) at a given 
value of Q 2 is determined by demanding minimal sensitivity of F n (M 2 , so) 
on the Borel parameter M 2 in the fiducial interval of the SR. These intervals 
M 2 G [M 2 /2, M?/2] are determined from the corresponding two-point NLC 
QCD SR and turn out to be M 2 = 1 GeV 2 , M 2 = 1.7 GeV 2 for the 
minimal NLC model and M 2 = 1 GeV 2 , M 2 = 1.9 GeV 2 for the improved 
NLC model, while the associated values of the pion decay constant read 
f n = 0.137 GeV 2 and f n = 0.142 GeV 2 , respectively. Notice that the value 
of the Borel parameter M 2 in the three-point SR roughly corresponds to the 
Borel parameter in the two-point SR, having, however, twice its magnitude. 
A stable window for the Borel parameter is obtained for thresholds in the 
range between 0.65 and 0.85 GeV 2 pQ. As a rule, the higher the value of 
sq, the larger the form factor because the perturbative input increases. The 
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sensitivity of the obtained form-factor predictions on the particular NLC 
model employed is rather weak in the region of momenta delimited in Fig. 
[T]by two vertical broken lines. We close this discussion by remarking that 
the overall agreement between the obtained predictions and the available 
experimental data from |18|, 1 19| is rather good. Moreover, they comply with 
the recent lattice calculation of [20], which is shown in the same figure in 
terms of a dark-grey strip, bounded from high-Q 2 values at Q 2 ~ 4 GeV 2 . 

In Fig. 2, we compare our predictions with other theoretical results, de- 
rived from the AdS/QCD correspondence (so-called holographic QCD), and 
some other models (including also the experimental data). As in the previ- 
ous figure, the shaded bands show our predictions for the minimal (dashed 
lines) and the improved NLC model (solid lines). The dashed (red) line at 
low-Q 2 (terminating at about 4 GeV 2 ) represents the prediction derived in 
[3] from QCD SR, while the thicker broken line below all other curves is the 
result of a calculation [2T] based on Local Duality QCD SR. The predictions 
from AdS/QCD are denoted by the upper short-dashed (green) line — soft- 
wall model [15] — and the penultimate dash-dot-dotted (green) line, obtained 
with a Hirn-Sanz-type holographic model [16] . The LD result of [21] is 
shown as a dash-dotted (red) line. 
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Fig. 2. Scaled pion form factor Q 2 F n (Q 2 ) within a band including uncertainties for the 
minimal (dashed lines) and the improved (solid lines) NLC model, using in both cases 
A 2 , = 0.4 GeV 2 . The experimental data are as in Fig. [T] The predictions from AdS/QCD 
are: upper short-dashed (green) line — soft-wall model ..15j dash-dot-dotted (green) line — 
Hirn-Sanz-type holographic model [IB] . The dash-dotted (red) line corresponds to the 
LD result of [21]. 

4. Conclusions 

Here we have studied a three-point AAV correlator within the QCD 
sum-rule approach with nonlocal condensates in order to obtain predictions 
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for the spacelike pion form factor pertaining to that momentum region ac- 
cessible to experiment at present and in the near future. The full-fledged 
analysis can be found in [1], where we also included in our discussion the 
so-called Local-Duality approach [5]. 

The principal ingredients of our approach are a spectral density that in- 
cludes 0(a s ) corrections, whereas the coupling used has an analytic struc- 
ture without Landau singularities, and an improved Gaussian ansatz for 
the nonlocal condensate. Our findings are supported by both the existing 
experimental data and also a recent lattice calculation in the momentum 
range up to approximately 10 GeV 2 . 
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